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Abstract 

Even in the absence of a sizable tensor contribution, a B-mode polarization can be gen- 
erated because of the competition between a pseudo-scalar background and pre-decoupling 
magnetic fields. By investigating the dispersion relations of a magnetoactive plasma supple- 
mented by a pseudo-scalar interaction, the total B-mode polarization is shown to depend not 
only upon the plasma and Larmor frequencies but also on the pseudo-scalar rotation rate. 
If the (angular) frequency channels of a given experiment are larger than the pseudo-scalar 
rotation rate, the only possible source of (frequency dependent) B-mode autocorrelations 
must be attributed to Faraday rotation. In the opposite case the pseudo-scalar contribution 
dominates and the total rate becomes, in practice, frequency-independent. The B-mode 
cross-correlations can be used, under certain conditions, to break the degeneracy by disen- 
tangling the two birefringent contributions. 



In the ACDM paradignj±| a potential candidate for the B-mode polarization are the ten- 
sor modes of the geometry inducing a frequency-independent polarization of the Cosmic Mi- 
crowave Background (CMB in what follows). By frequency-independent signal we mean that 
different observational channels measure angular power spectra with the same amplitude. 
In the opposite case the angular power spectra will effectively depend upon the (angular) 
frequency of observation. The WMAP 5-yr data pQ constrain the presence of a B-mode and, 
indirectly, r^, i.e. the ratio of the tensor power spectrum over the scalar power spectrum pQ. 
A further (frequency independent) source of B-mode polarization is cosmic shear (see e.g. 
[2]). Diverse data sets (such as the ones of Quad and Capmap [3]) impose concurrent limits 
on the B-mode polarization. Forthcoming experiments are expected to improve the present 
status of the observations by reaching into the region < 0.2. 

The only frequency-dependent signal investigated so far is provided by the Faraday effect 
which is a distinctive feature of magnetized plasmas in different contexts [4]. Large-scale 
magnetic fields present prior to the equality time are known to impact both on the temper- 
ature autocorrelations as well as on the polarization observables [5]. It has been recently 
shown, within a dedicated numerical approach [6J, that the Faraday rotation signal induced 
by a pre-decoupling magnetic field can overwhelm the B-mode polarization induced by the 
standard tensor contribution [7]. 

The B-mode autocorrelations might not be always sufficient to infer the presence of a 
pre-equality magnetic fielco. In short, the idea is the following. Consider a set-up where 
the pre-equality plasma is birefringent because of the concurrent presence of a pseudo-scalar 
background field, be it a, and of a large-scale magnetic field. Absent any pseudo-scalar 
background, the rotation rate would scale with the square of the wavelength [I] of the 
observational channel. In the presence of a pseudo-scalar field the dispersion relations can 
be generalized and the total rotation rate will have, both, a magnetic and a pseudo-scalar 
contribution. The purpose of this paper is to compute the B-mode polarization generated 
by the competition of the two aforementioned effects and to scrutinize if (and when) the two 
effects can be, at least partially, disentangled. The essentials of the problem at hand are 
usefully introduced in terms of the electromagnetic part of the action 



where g = detg^ and g^ u = a 2 (r)^; F^ u is the electromagnetic field strength; F^ u = 

1 Th.e A refers to the dark energy component (assumed to be in the form of a putative cosmological 
constant). The CDM refers to the (cold) dark matter contribution. In what follows the cosmological 
parameters will be fixed to the best fit of the WMAP-5yr data alone, i.e. (fibo> ^c0> ^A, ho, n s , e) = 
(0.0441, 0.214, 0.742, 0.719, 0.963, 0.087). In the latter string of parameters fl x denotes the critical fraction 
of a given species, ho fixes the present value of the Hubble rate; n s is the spectral index of curvature 
perturbations and e is the reionization optical depth. 

2 Following the established terminology the B-mode autocorrelations are denoted by BB. With similar 
notation we will talk about the TT, TE, EE angular power spectra meaning, respectively the autocorrelations 
of the temperature, the autocorrelations of the E-mode and their mutual cross-correlations. 
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e^ upa F pa / (2^—g) is the dual field strength in curved space-times. In Eq. ([T]) (3 is a coupling 
constant and M a typical mass scale which may take specific values, for instance, in a given 
scenario [5]. In Eq. ([I]) j u denotes the electromagnetic current which can be specified in 
terms of the charge carries (i.e. electrons and ions) as f = j\ +j% — e(n i -u 1 1/ — n c u^) (recall 
that, for both species, g^u^u^ = 1). The relevant set of equations can then be written, 
for brevity, in their covariant form and they are 

V M F^ = Anf + ^V^F^, = 0, (2) 

V M (7T + Tt u + T&) = 0, T% = p e>i < ; < ji5 (3) 

where is the energy-momentum tensor of the electromagnetic field and p ej i denote the 
energy density of electrons and ions. Note that V M (i.e. the covariant derivative associ- 
ated with the space-time geometry) does not only depend upon the scale factor but also 
upon the inhomogeneities. Equation ([3]) summarizes schematically the evolution equations 
of charged species whose governing equations can be more appropriately derived from the 
Vlasov-Landau equations in curved space (or from their lowest moments). While electrons 
and ions are coupled through Coulomb scattering, the electron-ion fluid is coupled to the pho- 
ton background. The plasma contains a large-scale magnetic field whose Fourier amplitudes 
satisfy 

(BiikrfBjfar)) = —P B (k)P t] (k)5^(k + p), P B = A B (—) (4) 

where k 2 Pij{k) = {k 2 5ij — kikj); A B is the amplitude of the magnetic power spectrum at 
the (comoving) magnetic pivot scale k^ (equal to 1 Mpc -1 in the forthcoming numerical 
examples). The magnetic field is inhomogeneous over typical length-scales which are of the 
order of the Hubble radius rn- The Larmor radius of the electrons, on the contrary, is 
tl — O^Vth/uBc) <^ r H where U Be is the (comoving) Larmor frequency and ~ \jT e /m e is 
the thermal velocity of the electrons. For a comoving field strength (9(nG) the Larmor radius 
is roughly eight orders of magnitude smaller than the Hubble radius. The guiding centre 
approximation (originally due to Alfven j9]) can be then applied. The charged particles 
orbiting around the magnetic field lines will see, in practice, a constant field up to drift 
corrections (going as [T C (B x V)B)/ (eB 3 ) where B = \B\) and curvature corrections (going 
as [T e Bx (B-V)B}/(eB 4 ) ) which are, however, negligible when the scale of (spatial) variation 
of the magnetic background is much larger than the gyration radius of the charge carriers. 
As discussed in [8] the dispersion relations can be derived by studying the propagation of the 
electromagnetic waves in a magnetoactive plasma at finite density. By writing Eqs. ([2]) (J3j) 
in their explicit form the compatibility of the system can be ensured if det^ljj = where Aij 
is given by 

= k 2 (Sij - kikj) - uJ 2 eij(uJ, a) + i-j^a 1 'e mij k m . (5) 

In Eq. ([5]) ejj (uJ, a) is the dielectric tensor and a' denotes a derivation with respect to the 
conformal time coordinate r. By requiring that det^j = and by setting the comoving 
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wavenumbers in such a way that k x = and k y = A; sin?? and k z = fccos?? the standard form 
of the Appleton-Hartree equation can be easily recovered and it is [8] 
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p <v 



where a = iTi/oJ, TC = a' /a and n = k/uJ is the refractive index. Denoting with uJ 
the comoving plasma frequencies for electrons and ions and with u;Be,i, the corresponding 
gyration frequencies e±(UJ,a) and e\\(uJ,a) are 
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The frequency U CT = /3a' /M measures the rate of variation of the polarization because of the 
presence of the pseudo-scalar background. If a is not homogeneous the dispersion relations 
will have a different form. According to Eq. ([6]) the refractive indices for electromagnetic 
propagation along the magnetic field (i.e. t? = 0) can be deduced from 
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(9) 



In the orthogonal direction the dispersion relations lead to the so-called ordinary and ex- 
traordinary plasma waves which are, however, not excited because of the minute values of 
the (comoving) plasma and Larmor frequencies for the electrons: 



u Be = 0.01759 
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(10) 



Indeed oJ max > a7 pe 3> aJse where aJ max = 27rz7 max and z7 max = 222.617 GHz corresponds to 
the maximum of the CMB spectral energy density. The rotation rate experienced by the 
linearly polarized CMB travelling parallel to the magnetic field direction is 
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Since uj a oc a 1 the contribution to the rate can be, in principle, either positive or nega- 
tive. This will affect the sign of the cross-correlations (e.g. the TB and EB angular power 
spectra). If to < uj a , the shift in the polarization plane of the CMB will essentially be 
independent upon the channel of observation §. In the opposite case (i.e. uj > u3 a ) the 
magnetized and the pseudo-scalar contribution concur in determining the total amount of 
rotation and, ultimately, the various polarization observables, i.e., according to Eq. ( TTTT) 
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Figure 1: The angular power spectra of the B-mode autocorrelations (plot at the left) and 
the absolute values of the cross-correlations (plot at the right) are reported in the case when 
the pseudo-scalar background and the magnetized background are simultaneously present. 
The ACDM parameters have been chosen in accordance with the best fit to the WMAP5-yr 
data alone pQ. 



ip{n : r) = < & ct (t) + ( &Faraday("-, T ) ■ The B-mode polarization induced by the magnetoactive 
plasma in the presence of a pseudo-scalar background can be computed by means of an itera- 
tive approach which generalizes the calculation of |7J. Since A±(n, r) = Aq(ti, t) ±z Au(nr), 
it transforms as a spin ±2 for rotations around a plane orthogonal to the direction of prop- 
agation of the radiation. The three-dimensional rotations and the rotations on the tangent 
plane of the sphere at a given point combine to give a 0(4) symmetry group [TU]. General- 
ized ladder operators raising (or lowering) the spin weight of a given function can then be 
defined as [IS EE]: 



K*Jh) = -(sintf) ±s 



sin it 



(sintf) TS , n = (#,¥>). (12) 



In real space the E-mode and the B-mode polarization will have spin weight s = 0: 

A E (n,r) = ~{K^(n)[K^(n)A + (n,r)} + Kt 1 \n)[K^ 2 \n)A.(n,r)]}, (13) 
A B (n,r) = l{K?\n)[K<?\h)A + (n,T)] - ^(n^K^^A^r)}}. (14) 

The heat transfer equation will then contain, in Fourier space, a convolution. Since the 
polarization is generated rather close to last scattering the iterative procedure of [7] implies 
that, to zeroth order in the rotation rate, the polarization is given, in real space, as 



A P (n,r) = ^3/2 J d 3 k Ap(fc, ji, r ), S P = A 



PO "I" Ap 2 + A12 



3 Different experiments are characterized by different channels of observations. For instance Quad [3] 
employs two series of bolometers located, respectively, at 100 GHz and at 150 GHz. 
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Ap(fc,^7- ) = 7(1 -/i 2 ) r JC(r)S P (k,T)e- ik ^-^dr, 
4 Jo 



(15) 



where /C(r) is the visibility function. In terms of Ap(n, r), Eqs. (fT3|) and (I14p imply 

A E (n, r) = -^{cos [2^(n, r)](l - /i 2 ) A P (n, r)}, 

A B (n, r) = 0j{sin [2^(n, r)](l - ^ 2 )A P (n, r)}, (16) 

where, as usual, /i = cos 1? and <9 M denotes a derivation with respect to cos In the absence of 
the (inhomogeneous) magnetized contribution Eq. ( jTBT) leads to the expressions customarily 
used in standard analyses, i.e. for instance 

C? B) = sin 2 2^f E) , Cf B) = -sin4^f E) , B) = sin2^f E) , (17) 

2 

^EE) (XE) 

where ip(n,r) = & a ( r ) — &a{ T o) is fully homogeneous and where the C e and C e are 
computed from Ap(n, r). If also the magnetized contribution is taken into account the 
situation changes both qualitatively and quantitatively. This aspect can be understood from 
Eq. (1161) . In the limit of small rotation rate 

A E (n, r) = - /i 2 )A P (n, r)], A B (n, r) = 2 0jty,(n, r)(l - ^)A P (n, r)}, (18) 

where ip(n,r) and Ap(n, r) depend upon the same point on the microwave sky. Eq. ([TBI 
holds in real space. In Fourier space the B-mode would be a convolution. In multipole space 
the angular power spectra inherit a peculiar form which contains a double sum involving 
also a known Wigner coefficient arising from the integral of three Legendre polynomials [7] . 

In Fig. [1] (plot at the left) the B-mode autocorrelations are reported in the case when the 
magnetized background competes with the pseudo-scalar background. The B-mode signal 
is frequency dependent. In the plot at the right the cross-correlations of the B-mode with 
the other CMB observables are illustrated. The observational frequency channel has been 
taken, for illustration, V = 100GHz. In Fig. CD denotes the magnetic field intensity 
regularized over a the pivot length-scale k^ 1 ~ Mpc. If the pseudo-scalar contribution is 
totally sub leading the Cf B angular power spectrum will diminish with the frequency as V~ A . 
Still the cross-correlations of the B-mode with the temperature and the E-mode polarization 
(i.e. Cj B and Cf B ) are non-vanishing. This aspect is illustrated in Fig. CE] (plot at the right) 
where Cj B and Cf B are reported. The numerical calculation leading to the results reported 
in Fig. [1] has been performed by including the magnetic field in the initial conditions and at 
every step of the Einstein-Boltzmann hierarchy and for the initial conditions corresponding 
to the magnetized adiabatic mode. 

In summary, if observations point towards a frequency dependence of the B-mode po- 
larization Faraday rotation is probably the only candidate. The cross-correlations (i.e. the 
EB and TB spectra) are expected to vanish in the case of a stochastic magnetic field leav- 
ing unbroken spatial isotropy [5]. If they are observed this means that the rotation rate is 
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quasi-homogeneous and a pseudo-scalar background field may be around. In the latter case, 
if uJ > uJcr the B-mode will scale with frequency as dictated by the Faraday effect while the 
EB and TB correlations will allow to measure independently u a . In the opposite case (i.e. 
ZJ < uJa) the frequency dependence induced by the Faraday effect is overwhelmed by the 
(homogeneous) pseudo-scalar rate. In this second case the effects of the primordial magnetic 
fields will be imprinted on the EE and TT angular power spectra [5J E] but the B-mode 
autocorrelations will be independent of the frequency. This demonstrate that the analysis 
of the B-mode autocorrelation is necessary but might insufficient, if taken individually, to 
infer the existence of pre-decoupling magnetic fields. 
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FIS2006-05319 and CSD2007-00042 of the Spanish Science Ministry. 
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